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Abstract 

It is discussed the possibility of a fine-tuneable contribution to the 
two way Doppler acceleration either towards, either outwards the Sun 
for heliocentric distances above 20 AU by considering a background 
described by an Expanding Locally Anisotropic (ELA) metric. This 
metric encodes both the standard local Schwarzschild gravitational ef- 
fects and the cosmological Universe expansion effects allowing simulta- 
neously to fine-tune other gravitational effects at intermediate scales, 
which may be tentatively interpreted as a covariant parameterization 
of either cold dark matter either gravitational interaction corrections. 
Are derived bounds for the ELA metric functional parameter by con- 
sidering the bounds on the deviation from standard General Relativity 
imposed by the current updated limits for the Pioneer anomaly, taking 
in consideration both the natural outgassing and on-board radiation 
pressure, resulting in an average Doppler acceleration outwards the 
Sun of Op w +^-^-2.0 ^ 10"-^° (ms~^). It is also computed the mass- 
energy density for the ELA metric within the bounds obtained and are 
discussed the respective contributions to the cosmological mass-energy 
density which, for compatibility with the AC DM model, are included 
in OcDAf. 



keywords: gravitation, Pioneer, interplanetary medium, equation of state, 
dark matter, cosmological parameters 



1 Introduction 



The Pioneer space-crafts have been launch from Earth on 1972 and 1973 
travelhng outwards of the Solar system [TH3]. In between the heliocentric 
distances of 20 AU and 70 AU, travelling at an approximately constant ve- 
locity of Vp = 12.2 X 10'^ ms^^, it was detected an unmodeled Doppler shift 
known as the Pioneer effect. Originally attributed to a physical acceleration, 
this effect was interpreted as a constant radial acceleration towards the Sun 
of Op = —(8.74 lb 1.33) X 10^^^ m [3]. Many possible explanations for 
this effect have been considered in the literature, including modified gravity 
theories and models |4j, non-gravitational interactions [5j, extra-dimensional 
theories and models [6] and dark matter [?]• 

However recently the Pioneer data was re-analysed favouring a significant 
jerk in the measured Doppler acceleration such that the best exponential fit 
is approximately given as a function of the radial distance to the Sun (for 
r e]20 AU, 70 AU[) by 

.n-ir- ir, 2.599x10^-'-^ log(2) / , .ncs,inll\ r, 

ae.p{r) ~ -12.22l0:|6 x IQ-^^e ranj^ ('-i-496xio") ^^^-2) ^ 

In addition it has been put forward that this acceleration is mainly due 
to two distinct effects. The more significant contribution ay is of thermal 
origin, the antenna located at the back of the space-craft works as a radiation 
sail propelled by the radiation emissions from the on-board power source 
and the space-craft electronics [TUHT^ . Also the natural outgassing of the 
space-craft material [12] contributes a significant amount ai for the Pioneer 
acceleration. These two contributions are (for r €.]20 AU, 70 AU[) |12j 

arir) « -10.7t|l X 10-10e-6-238xl0-^(r-1.496xl0") (^^-2) ^ 

(2) 

aUr) « -7.15tlfo x lO'ii (l + 2-389x10^ ^ (ms'^). 

The resulting acceleration in the range r g]20AU, 70AU[ is, within the ex- 
perimental plus the estimative error bars, null. Its average value in this 
range is 

(ap(r)) = (ae,p(r) - aUr) - aT(r)) « 0.4^|i x 10"!° (m^-^) . (3) 

The several contributions for the unmodeled acceleration of Pioneer 10 and 
the resulting acceleration are plotted in figure [TJ 

The main purpose of the this paper is to investigate whether a contri- 
bution to the Doppler acceleration of the Pioneer space-craft can be con- 
sistently obtained considering a background described by the ELA met- 
ric an anisotropic [TJ] ansatz that generalizes the isotropic McVit- 
tie metric |16H19] interpolating between the Schwarzschild (SC) metric [20] 
near massive bodies and the Robertson- Walker (RW) metric ^21j at spatial 
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(a) Contributions to the acceleration of Pioneer 10 



(b) Excess Acceleration of Pioneer 10 
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Figure 1: (a) contributions aL[r) and axir) ^ to the unnnodeled acceleration 
<iexp{f) (QP of Pioneer 10; (b) total Pioneer 10 acceleration ap{r) = aexp{f) — 
aL{r) — aT{r). The shaded region represents the quoted experimental plus estimative 
error bars JMESMM- 



infinity, hence describing local matter in an expanding Universe [22]. In 
particular if the bounds on the Pioneer 10 acceleration discussed above can 
be employed to constraint the parameters of such metric. 

The shift function of the ELA metric can be fine-tuned through an ex- 
ponential functional parameter allowing for a covariant parameterization of 
intermediate scale gravitational effects |14ll23j . hence allowing to describe 
either local cold dark matter effects |24y25j either modified theories of grav- 
ity [26] . Although this fine-tuning is technically more challenging than semi- 
classical heuristic approaches it has the following advantages: 

• it is a covariant formulation allowing to deal with gravitational cor- 
rections within the framework of General Relativity; 

• in particular allows for direct covariant computations of gravitational 
accelerations and gravitational frequency-shifts; 

• allows as well for QFT cross-section computations in curved back- 
grounds to include cold dark matter effects; 

• it allows for the description of intermediate scale gravitational effects 
maintaining the local and cosmological physical properties of space- 
time. 

In the following we will employ two distinct coordinate systems, expand- 
ing coordinates, ri, for which the area of the sphere is ^ = Anrf such that 
the coordinate quantities correspond directly to measurable lengths and the 
standard Robertson- Walker coordinates, r, for which the area of the sphere 
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is A = Att a^r'^ being suitable for direct computations of measurable red- 
shifts [U]. The map between both these coordinate systems is ri = ar and 
the ansatz for the ELA metric is |13p i4] 

= {1 - U){cdtf - (inl 



(4) 

= {l-U){cdtf -a^d^'^ 

where a is the universe scale factor, H = a/a \s the time dependent Hubble 
rate, U = 2GM/{c^ri) = 2GM/{c^ar) is the standard SC gravitational 
potential, M is the standard SC mass (the value of the gravitational mass 
pole), c is the speed of light, G is the Gravitational constant and the solid 
angle line-elements are dQ^ = rf{d9'^ + sir? 6 d(p^) and dVl^ = r'^{d9'^ + 
sin2 edf^). 

This ansatz generalizes previous ansatz and solutions in the literature [161 - 
[19] allowing for local anisotropy [15] and having the novelty of describing a 
space-time without singularities at the SC radius simultaneously maintain- 
ing as asymptotic limits the SC metric near the SC event horizon and the 
RW metric at spatial infinity. To ensure these properties of space-time the 
following limits must be obeyed [TBlfll] 

lim (1 — = : SC mass pole M is maintained , 

ri--s>0 



lim a > 3 : space — time is singularity free 

at SC horizon 

lim (1 — ?7)" = 1 : metric converges to RW metric 



(5) 



The main drawback of this metric is the absence of a direct physical 
interpretation for the functional parameter a. Here we will not attempt to 
solve this problem, we recall that the ELA metric was originally introduced 
phenomenologicaly in |131ll4j to generalize the McVittie metric [16J and it 
is at most interpreted as a heuristic fine-tuneable parametrization of grav- 
itational interactions at intermediate scales. Hence we consider a branch 
ansatz with the following properties: 

• allows the limits (0) to be obeyed both near the origin and at spatial 
infinity; 

• has no relevant contributions for planetary motion within the solar 
system; 
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• allows to fit the Doppler acceleration for the heliocentric distance of 
20 AU up to 70 AU; 



• allows to model the mass-energy density described by this background 
for distances above 70 AU up to some upper cut-off distance Rmax- 

In particular, above the cut-off Rmax-, we explicitly consider the McVittie 
metric, hence a = 0, such that explicitly isotropy is retrieved for large radial 
distances and the mass-energy density for the background metric coincides 
with the one for the RW metric. Also we do not discuss here possible effects 
on planetary motion, leaving an account of this problem to a future work. 
Here for heliocentric distances below 20 AU we consider the lower possible 
value of a ~ 3 compatible with the limits ([5]) such that the gravitational 
effects due to the background are negligible for planetary dynamics within 
the inner Solar system. Therefore we consider the following expression for 
the functional parameter a: 



where — l^Qi<Oisa negative constant arbitrarily close to zero that 
ensures that the SC mass pole is maintained at the origin [13] being its effects 
negligible far from the SC horizon, a{2GM/c^ < ri < 20 AU) w 3. As for 
02 (r) is a generic function that can be fitted to the Pioneer experimental 
data above the heliocentric distance of 20 AU . 

We have organized this work as follows. In section [2] are derived the 
Doppler shift for a two way range Doppler measurement as well as the grav- 
itational acceleration correction for the background described by the ELA 
metric ^ . Then in section [3] are derived the mass-energy density and the 
pressures as well as the mass corrections due to the extended configurations 
corresponding to the functional parameter ansatz It is also discussed 
the respective contributions to the relative cosmological mass-energy densi- 
ties within the AC DM model [27]. In section d] are computed and discussed 
several fits to the functional parameter a within the error bars of the Pio- 
neer 10 data. Finally in the conclusions are resumed the results obtained 
and discussed its relevance as well as future research topics. 

2 Doppler shift and gravitational acceleration 

The Pioneer effect was measured using a two-way Doppler shift [HIS]. The 
standard gravitational red-shift as described by the SC metric is null, as 
for the standard cosmological red-shift is negligible as noted in [3l[Tl] being 
below the accuracy of the Pioneer measurements. For this reason the Pioneer 



' {3-ai) + aiU 





max 



(6) 
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effect is often attributed to a correction to tlie gravitational acceleration [U 

ElEj. 

Taking as central mass the Sun and neglecting the contribution of other 
bodies in the solar system, the corrections to the gravitational acceleration 
felt by a test mass travelling in the background given by the ELA metric dH 
with respect to the General Relativity gravitational acceleration correspond- 
ing to a Ricci-flat background described by the SC metric is 



K « ^ (l-C/0r(2(l-C/e)-(l + a)r,C/^ 



l_a 
I 2 2 



-2{1 + q){l - Uq)- 
+rp(l - Uq) ln(l - UQ)a') + 0{H' 



(7) 



where the primed quantities stand for differentiation with respect to the 
radial coordinate ri, stands for the geometrical (measurable) distance 
from the Sun to the Pioneer space-craft and Uq = 2GM(7)/{c^rp) stands for 
the standard SC gravitational potential of the Sun with Mq = 1.9891 x 
10^° kg. 

In addition there is a correction to the gravitational red-shift. The mea- 
surable red-shift due to the background is obtained from the line element (jH) . 
For a light ray travelling approximately at a radial trajectory with respect 
to the Sun described by the four- vector k^^ k. (uj/c, k"^ ,0,0) , from earth at 
r^f^ = gqA^^ ~ 1 AU to the Pioneer space-craft at r^^ = ai r^^'^ = rp and 

back to earth at r^^ = 02 r*-^) ~ 1 AU we obtain that the ratios between the 
radial wave numbers of the emitted radiation from earth k'^^^ , the radiation 
received at the Pioneer space-craft k^-^-^ and the radiation received back at 

^[1) _ ao 1-Up- 5 



Earth kj^^-^ are 



e 



^(2) ^ Ql 1 - (8) 

a2 I -Up + 5 

S = _[/)§+ 

c V 

where Up and Ue are the gravitational potential of the Sun evaluated at the 
Pioneer space-craft and at Earth, and gq, ai and 02 stand for the Universe 
scale factor a = a{t) evaluated at the time of emission of the radiation 
from Earth, the time of receiving and re-emission at the Pioneer space-craft 
and the time of receiving it back at Earth, respectively. Further noting 
that expanding the Universe scale factor to first order in the Hubble rate we 
obtain that oi = ao{l+iH/c)rp)+0{H'^) and 02 = ao{l+2{H/c)rp)+0{H'^), 
the measurable frequency-shift (red-shift if positive or blue-shift if negative) 
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due to the background for a two way Doppler measurement from Earth to 
Pioneer and back is 




1-Up-S 



1 + 0{H^) 



(9) 



1 + 2 



lllR l-U„ + 6 



Hence the contribution to the Pioneer non-physical acceleration due to the 
background red-shift is obtained by differentiating the standard two way 
Doppler shift Au/vq = —2fp/c 



where the prime denotes differentiation with respect to the radial geometric 
distance (the coordinate ri), Vp = rp for the radial geometric velocity of 
the space-craft (obtained from the differentiation dt = dtvi dr^ ) and we have 
neglected relativistic corrections (7 ~ 1 for <C c). 

Before proceeding to fit these results to the Pioneer acceleration bounds 
represented in figure [1] we will further derive and discuss the mass-energy 
density of the background described by the ELA metric (jH). 




(10) 
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3 Mass-energy density and pressures 



The mass-energy density, and anisotropic pressures Pr(a) ^-nd Pe(a) 
Ptp{a) ^or the ELA metric ^ are 



Pia) = ^{l-Ur-'{aU+{l-U)x 



X (3 + rilog(l-C/)a,,a)) , 



Pria) = ^{l-U)"^ (2{l + q)-{l-U)'^{aU 



+ (l-C/)(3 + rilog(l-f/)a,,a) , 



4(l + g)(l-C/(l-^)) (11) 



.2a-.)t^(3.>-f.(i-|))(i-|) 

-3 - riC/(l - U)dr,a) + ri(l - ^7) log(l - ^7) 



x{^{l + q)dr^a -2U{1- U)~ adr^a 
-{l-U)"^ {Qdr.a + n log(l - U){dr,af 

When the hmits jS]) are obeyed [13], at ri = ri.sci P(a) coincides with the 
mass-energy density for Ricci flat space-times, being null — > ?^i.sc) = 
0, and at spatial infinity it asymptotically converges to the RW mass-energy 
density — )• -|-oo) = pRw = 3/i^/(87rG), as well as +oo) = 

Pe{a)iri +oo) = p^(a)(ri +oo) = pRw = (2g - 1) /i^/(87r G) such that 
at spatial infinity the RW equation of state is obtained ojrw = (2^ — l)/3 = 
—0.728 (this value is obtained within the ACDM model assuming that the 
pressure of the background is given only by the cosmological constant [27]). 
Hence this mass-energy distribution has been interpreted in [I3J , for a con- 
stant functional exponent a as a local anisotropic extended correction to the 
global isotropic cosmological background due to a local point mass vanishing 
at spatial infinity. 

In addition to preserve causality it is required that the mass energy 
density be positive, hence /9(a) > 0. This inequality is bounded by the 
solution of the differential equation a{ri) = —(1 — [/(ri))(3 -|- rilog(l — 
U {ri))a' {ri)) /U {ri) holding the following inequality 

^ ^ aconst - 3 log (n) 

"^^^^ - log(i - U[n)) ' ^''^ 
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where O-const is an arbitrary integration constant. 

With respect to the contribution to the cosmological mass-energy den- 
sity due to the ELA metric, it can be estimated by comparing the total mass 
contribution of the mass-energy density /0(q) (jlip with the central mass con- 
tribution Mq. To compute such contribution it is enough to integrate the 
difference between the and the RW expanding background mass-energy 
density pRw = 3 x H'^ /{SirG). Specifically for the branch ansatz d^D we 
obtain 

r+oo 

Ma = 47r / rl{p(^a) - pRw)dn 
Jo 

+ ((l-?7(20AC/))"(2°~^^)-l) {20 AUf 



If ^ 
2G L 



(13) 



1 - C/(20AC/))"(2°^^^) - l) {20 AU) 

where the last two lines are the contributions due to the discontinuity at ri = 
20 AU and above R 

max the ansatz exactly matches the McVittie metric (a — 
0) such that = pRw and the mass contribution is null. We note that as 
long as the limits ([5]) are obeyed near the origin, besides the standard SC 
mass pole Mq , there is no further mass contribution due to the divergence at 
the origin. In addition convergence of Mq, when the cutoff Rmax is taken to 
infinity {Rmax — > +oo) requires more severe conditions than the limits ([S]), 
in particular it can be fine-tuned either to be null (e.g. for a2{ri — )■ -|-oo) = 
constant/rj^, n > 2) or to a multiple of the central mass (e.g. for 02(^1 — ?■ 
+CX)) = constant/rf we obtain Mq, = i?^ constant Mq/(2G)). Nevertheless 
keeping Rmax finite in the ansatz ^ keeps the quantity Mq finite and allows 
to fine-tune the total mass of the background. 

For exemplification purposes the contribution to the cosmological mass- 
energy density due to Mq, can be estimated by assuming that all barionic 
matter has a proportional contribution to the one obtained here for the Sun 
such that r^Q/rif, = Mq/M0, hence 

^^a = ^f^6, (14) 

where $7?, = 0.0456 is the relative cosmological barionic mass-energy within 
the AC DM model [27]. The contributions to Of, are given mostly by de- 
tectable massive bodies in the universe, hence a direct estimation is possi- 
ble. Also we recall that within the ACDM model the remaining contribu- 
tions to the total mass-energy density are due to the cosmological constant 
flA = —{2q — l)/3 = 0.728 and due to non detectable cold dark matter 
^CDM = 0.226 (also denominated J7c)- The radiation density as well as 
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the curvature density are negligible such that the total relative density is 
fl = Qf, + VtcDM + ^^A = 1- With respect to the cosmological equation 
of state we note that both the barionic matter and cold dark matter have 
null pressures, while the cosmological constant has a negative pressure with 
equation of state loj^ = — 1. Then the only contribution to the cosmological 
pressure is due to the cosmological constant such that the equation of state 
of the Universe \s u = pj^/Q. = —57 a = (2(? — l)/3 = —0.728. Further noting 
that, due to the upper cut-off Rmax in the ansatz for a ([6]), the average con- 
tributions to the cosmological pressures due to the ELA metric background 
are null 



pRmax 

Att / pr dr 

rRmax 

4tt / pg dr 



(15) 



we conclude that, within the ACDM model, contributions to the cosmolog- 
ical densities from the background @ with the functional parameter anzatz 
given in ([6]) must me included in ^Icdm- 

Possible alternatives to the ansatz that consider non-null values (a / 
0) for the functional parameter up to spatial infinity would also allow for 
non-null contributions to the cosmological pressure hence also allowing for 
a contribution to both the cosmological constant density and pressure. We 
do not develop this construction here. 

As for the local equations of state for the ELA metric we note that 
generally the pressures Pr{a) a-nd pe(a) = Ptp{a) a-^e negative far from the SC 
event horizon, in particular for constant values of the parameter a > 3, there 
are regions where both iOr[a) = Pr{a)/ P{a)) and oJei^a) = Pe(a)/{c^ P{a)) are 
above —1 and regions where u^^a) is above —1 and oj0(^ci) = Pe(a)/ {c^ P(a)) is 
below — 1 [13] . The interpretation of the matter content for these two regions 
is distinct, in the first regions the matter distribution may be interpreted as 
due to a local configuration of a standard scalar field (wr'(a) > — 1 and ujQ(^a) = 
^tf(a) > ~1) pS], while in the second regions the matter distribution is 
compatible with a standard scalar field along the radial coordinate direction 
(ujr > — 1), along the angular variables directions it is obtained an equation 
of state commonly attributed to Phantom matter (w < —1) [29]. We recall 
that the main difference between these two scalar fields is the sign of the 
kinetic term in the action, hence a implementation for these regions requires 
that the kinetic term along the angular variables directions to be negative 
while being positive along the radial coordinate directions. Possibly the 
simpler implementation for these equations of state is to consider a local 
configuration for a radially symmetric anisotropic scalar field. Such model 
may be compatible, for instance, with locally anisotropic string theories and 
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string inspired cosmological models [30j . 

Next we will employ the Pioneer bounds on rp to obtain bounds on 
the parameter a above the heliocentric distance of 20 AU and discuss the 
total contributions of such configurations to the cosmological mass-energy 
densities. 

4 Bounds on a 

A contribution to the Pioneer Doppler acceleration due to the background 
encoded in the metric dU is the sum of the contribution ([7]) and ([TO]) 

a, = 6rp + Srl^°fP''' . (16) 

Next we consider and analyse 4 distinct fits of the functional parameter 
o:2{i"p) within the acceleration bounds for ap as represented by the shaded 
region in figure [D 

o^max the allowed maximum valued fit for a corresponding to positive Doppler 
acceleration aQ,^^^ > as close as possible from the upper boundary 
of the allowed region; 

a+ a fit for a corresponding to positive Doppler acceleration > half 
way the positive acceleration allowed region; 

a- a fit for a corresponding to negative Doppler acceleration aa_ < half 
way the negative acceleration allowed region; 

amin the allowed minimum valued fit for a corresponding to negative Doppler 
acceleration aa^.^ < as close as possible from the from the lower 
boundary of the allowed region. 

First let us note that for positive values of the function 02 (r^) the effects 
are of the same order of magnitude of the pure background expansion, Cq, ~ 
+10~^^ms~^, hence to obtain measurable effects it is required to consider 
negative values for this parameter. 

We remark that an exact fit to the upper boundary of the allowed pos- 
itive Doppler acceleration Cq, correspond to strictly negative values for the 
mass-energy density /9(q,) in the range Vp g]20 At/, 70 AC/[, therefore for the 
fit (y-max < we consider the limiting solution (fT2|) with aa^^^{2QAU) = 
3.057 X 10^^*^7715^^, hence with oiconst = 106.176. As for the fit Cq,^, for 
exemplification purposes and to ensure positive p(^a) we consider the fit- 
ting function to be 9/10 of the limiting solution ([12]) with aa^^^{2{)AU) = 
1.529 X 10~^''ms~^, hence with aconst = 107.650. As for the a_ and amin 
corresponding to negative values for the Doppler acceleration Oq, < we 
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consider a power series fit given by the following expression 

N 

a2{rp) = ^ a2.nrp . (17) 

n=-l 

For the specific case of the a_ and Omm we consider A'" = 10 and the specific 
coefficients are listed in appendix [XI The 4 fits considered are pictured in 
figure [2] (a) . Noting that amin and Umax define the lower (negative) and 
upper (positive) limits on the Doppler acceleration allowed by the bounds 
on the Pioneer 10 acceleration represented in figure[T]we conclude that in the 
range rp g]20 AU, 70 AU[ the functional parameter 02 has the lower bound 
represented in figure [2] (b). 



(a) Allowed fitting region (b) Allowed values for itj 




20 30 40 50 60 70 20 30 40 50 60 

t, CAU) r, (AUl 



Figure 2: Fit of the functional parameter a2{ri) to the Pioneer accelera- 
tion bounds of figure (a) amax corresponds to the limiting solution with 
oiconst — 106.176, the solution a+ to 9/10 's of the limiting solution with 
oiconst — 107.650, the solution and amin correspond to a power series fit of 
the parameter a |J7[ ); (b) lower bounds for a < 0. 



The dominant contribution to for the above fits between heliocentric 
distances of rp g]20 AU, 70 AU[ is the background frequency-shift, 5rp°^^^^^ r-. 
_2Q-iOj^g-2^ As for the value of the physical gravitational acceleration, it 
is 4 to 5 orders of magnitude below this value, for Umax and a+ we obtain 
(Jr'p ~ — 10~^^ ms~'^ (towards the Sun) and for amin and a_ we obtain Jr^ ~ 
+10""*^^ ms~'^ (outwards the Sun). These values are well below the standard 
General Relativity acceleration which is of order vqh ~ —10~® ms~^. Nev- 
ertheless the effects of this correction to the gravitational acceleration can be 
measurable at least for Pluto and Neptune [31] , for instance if this correction 
is fully accounted by an orbital radius variation we would obtain at an orbital 
radius of Vorb = 20 AU for amax and a+ the value rorb/forb ^ — 10~^^ s^^ and 
for amin and a_ the value ri/vi ~ -|-10~^^ s~^. These values are well below 
the experimentally measured orbital radius variation for inner planets in the 



11 



Solar system such as Venus and Mars, \'rorb/i^orb\ < 10^^^ s^^ [32]. We post- 
pone a detailed account of the effect of the ELA metric in planetary orbits 
to another work, here we have considered that below the heliocentric dis- 
tances of 20 AU a(ri < 20 AU) = 3 ([6]) for which the physical gravitational 
acceleration is of order Vp ~ 10~^^ms~^ such that rorb/forb ~ — lO^^'^s^^, 
hence of the same order of magnitude of pure expansion effects. 

As for the mass-energy densities we note that up to the heliocen- 
tric distance of 20 AU it is slightly below the cosmological quantity prwi 
hence contributing a negative amount to the cosmological mass-energy den- 
sity (fH|) . /f3(Q,) — pbw < 0. Above the heliocentric distances of 70 AU the 
mass-energy densities for the fit o-max is exactly null due to saturating the 
bound (fT2]) while the remaining 3 fits, «+, a_ and amin correspond to 
mass-energy densities well above the cosmological quantity pi^a) ~ PRW ^ 0, 
hence while Omax contributes a negative amount to the cosmological mass- 
energy density, a+, a_ and amin contribute a positive amount. Although 
the Pioneer experimental data only extends up to the heliocentric distance 
of 70 AU, for exemplification purposes of estimating maximal contributions 
to the relative cosmological mass-energy density SIq, we extend the fits ob- 
tained up to some heliocentric distance Rrnax ^ 70 AU . Above Rmax 
consider a = (0) for which the ELA metric coincides with the isotropic 
McVittie metric [16] such that spatial isotropy is exactly recovered and 
> Rmax) = PRWi hence the contribution to Ma above Rmax due 
to the ELA metric background is null. The contributions for the 4 fits of 
the functional parameter a represented in figure [2l as well as the respective 
mass ratios Ma/MQ, are plotted in figure [3] up to the heliocentric distance 
100 AU. 




Figure 3: (a) mass-energy densities p^^) for the fits amax, ol+, a_ and amin 
represented in figure\^ (b) respective mass ratios M^/Mq, the negative mass ratio 
Ma^^^/MQ ~ — 10~^^ is several orders of magnitude below the remaining fits not 
being represented in the figure. 
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Therefore, in the range of the hehocentric distances rp €.]20 AU, 70 AU[, 
the fits obtained contribute at most an amount of order ~ 10~^ to the rela- 
tive cosmological mass-energy density Q,. We note that, as aheady discussed 
in the previous section, within the ACDM model these contributions can 
only be included in the relative cosmological mass-energy density of cold 
dark matter 0,cdm = 0.226. Hence these contributions are negligible up to 
70 AU. 

For exemplification and analysis purposes, let us extend these fits above 
the Pioneer data range of 70 AU. For the fit amax the mass-energy density is 
null, hence the contribution to is negative being given by —^b^§^^^ ^max (O- 
This contribution can be fine-tuneable being unbounded from below, hence 
increasing the required value of contributions to ^cdm of other sources 
that not the ELA metric background. For the fit a_ the mass-energy 
density is positive above the heliocentric distance of 20 AU up to spatial 
infinity decreasing monotonically with rp such that at rp = 10930 AC/ it 
exactly matches the cosmological quantity prw and vanishes at spatial in- 
finity. Hence the maximum positive relative contribution to the cosmological 
mass-energy density is Qa = +6.412 x 10~^^Qf) corresponding to the cut- 
off Rmax = 10930 AU. As for the maximum negative relative contribution 
to the cosmological density it is fine-tuneable and unbounded from below 
similarly to the fit amax- As for the fit a_, the mass-energy density pa_ 
is positive above the heliocentric distance of 20 AU up to approximately 
ri = 95.48 AC/ becoming negative above these distances such that the max- 
imum allowed value of the ansatz cut-off ([6]) coincides approximately with 
this value, Rmax ~ 95.48 AC/. For this fit the contribution to the rela- 
tive cosmological density is always positive having the maximum value of 

= +3.367 X 10-^°Ob for Rmax = 95.477 AC/. Finally for the fit Omin 
the mass-energy density Pa„^i„ is positive above 20 AU increasing monotoni- 
cally above the heliocentric distances of 96.596 AC/. Hence the contribution 
to the relative cosmological mass-energy density is unbounded from above 
exactly matching the cold dark matter contribution = 50;, = ^cdm for 
Rmax = 183.328 AC/. These results are summarized in table [TJ 

From the above discussion we conclude tlictt the fits ttmaa; 

and a+ rep- 
resented in figure O corresponding to a positive Doppler shift in the range 
rp g] 20 AC/, 70 AC/ [, are compatible with the cosmological ACDM model 
only when the fits are extended up to ~ 10 pc. When higher ranges for 
these fits are considered the negative contributions to the relative cosmolog- 
ical mass-energy densities become relevant such that the required relative 
(positive) density of cold dark matter to offset these contribution increases 
significantly hence becoming non-compatible with estimative for dark matter 
density in the universe [MIEZ]- As for the fits a_ and Umin, corresponding 
to a negative Doppler shift in the range rp g] 20 AC/, 70 AC/ [ have more re- 
strictive bounds being valid only up to order ~ 10^ AU. Above this range 
the mass-energy density corresponding to the fit a_ becomes strictly nega- 
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= 


fine-tunable 


— Qu^^f'^'^ R'^ < 

— '5 '6 3 jv/q ^max ^ 




> Prw 

(< PRW 


= 10930 AC/ 
> 10930 AC/ 


= +6.412 X 10-^°f]b 
< +6.412 X l{)-^^Q.b) 


a- 


> PRW 


= 95.477 AC/ 


= +3.367 X 10-iOf^b 




> PRW 


= 183.328 AC/ 


= +50fc = VtcDM 



Table 1: Resume of maximal contributions to the relative cosmological mass- 
energy densities for the fits to Umax-, «+> «- and amin- 



five and for the fit amin the relative mass-energy contribution exceeds the 
one attributed to cold dark matter. Hence the fits amax and 0+ allow for a 
relatively larg er range than the fits o^— and o^mm' 

For last let us discuss the equation of state corresponding to the func- 
tional parameter fits. We have already concluded that the cosmological 
equation of state is not affected by the ansatz for a dG]), the average correc- 
tion to the pressure is null (|15p and within the AC DM the corrections to the 
mass-energy may be included in the relative density of i^cDM- However lo- 
cally there are significant deviations from the cosmological equation of state 
as already analysed in [13] for a constant parameter a. Specifically for the 
limiting case corresponding to the fit amax the mass-energy density is exactly 
null from heliocentric distances of 20 AC/ up to Rmax, hence it is not possible 
to define a equation of state. As for the anisotropic equations of state cOr and 
Ljg = uj^p for the fits a+, a_ and amin are plotted in figure HI figure [5] and fig- 
urelH respectively. For the fit a+ the radial equation of state Ur above helio- 
centric distance of 20 AC/ is above —1, being negative up to 12154.3 AC/ and 
positive above this distance while the angular equation of state ojq = co^ is 
always positive. For the fit a_ the radial equation of state Ur is negative and 
above —1 in the range rp g]20 AC/, Rmax = 95.48[ while the angular equation 
of state is positive in this range. For the fit amin the radial equation of state 
ojr is negative and above — 1 in the range rp (z\'20 AU, Rmax ~ 

183.328 AC/[ 

while the angular equation of state is negative being below —1 in the range 
rp g]20 AC/, 25.86 AC/[U]96.60 AC/, i?,„ax = 183.328 AC/[ and above -1 in the 
range rp g]25.86 AC/, 96.60 AC/[. 

As for the possible interpretation for the matter content corresponding 
to the several fits discussed we note that the fits and a_ have equation 
of state compatible with a standard gauge field {tdr(a) > — 1 and tde{a) = 
<^</3(a) > ~1) while the fit amin has a equation of state along the radial 
direction compatible with a standard gauge field {ujj.(a) > —1) as well as 
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(a) U), for the fit «■+ 



(b) oif,=(jy^ for the fit n+ 



-0.9990 - 
-0.9992 - 
-0.9994 - 
-0.9996 - 



3* 0.3503 - 



20 30 40 50 60 70 80 90 100 
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Figure 4: Equation of state for the fit a+ of the functional parameter a2(^i) (0)' 
(a) radial equation of state uJr,' (h) angular equation of state ojq = oo^. 



(a) bj, for the fit a. 



(b)(u„=iu^for thefitit_ 




r,, (AU) 



Figure 5: Equation of state for the fit a- of the functional parameter a2(^'i) (0)' 
(a) radial equation of state uJr,' (h) angular equation of state ojq — oj^. 



along the angular variable directions in the range Vp g]25. 86 ^C/, 96.60 
while in the range g]20^[7, 25.86 AC/[U]96.60 At/, = 183.328 AC/[ it 

has a equation of state that resembles Phantom matter = uj^^^^-) < 

-!)• 



5 Conclusions 

In this work we have computed the measurable Doppler acceleration cor- 
rections for a background described by the ELA metric (H} and fitted an 
ansatz for the functional parameter of this metric ([6]) to the current Pio- 
neer anomaly bounds on the range Vp g]20^C/, 70 AC/[, hence setting lower 
bounds for the functional parameter a in this range. 
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(b) (uh=(u^ for the fit 



20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 1 
r,, (AUl 




20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 ISO 
r,, lAU) 



Figure 6: Equation of state for the fit amin of the functional parameter a2{ri) 
(a) radial equation of state iOj.; (h) angular equation of state ujg — u)^. 



In particular we have consider the 4 fits represented in figure [5] within 
the error bars for the Pioneer Doppler acceleration, the fit Umax correspond- 
ing to the limiting solution (jl2p . the fit a+ corresponding to 9/10 of the 
limiting solution ()12p and the power series (|17p fits a„ and amin- The 
main contribution to the Doppler acceleration corresponding to these fits 
is due to the background gravitational frequency-shift (jlOp . being of order 



6r. 



Doppler 



10 



-10 



ms 



^ while the physical gravitational acceleration ([7]) is 



lower by 4 orders of magnitude, 5fp 

We have also computed the mass-energy density and anisotropic pres- 
sures for the ELA metric (jlip and analysed the possible contributions due 
to this background to the relative cosmological mass-energy densities and 
pressures. There are no contributions for the cosmological pressure, however 
extrapolating the fits above the heliocentric distance of 70 AU are obtained 
significant contributions to the cosmological mass-energy density. For com- 
patibility with the AC DM model we conclude that these corrections must 
be included in the mass-energy density attributed to cold dark matter Qcdm 
such that the cosmological equation of state is not modified. Hence we note 
that the ELA metric background may constitute a covariant parameteriza- 
tion for cold dark matter. We also discuss the anisotropic local equation 
of state for the fits obtained noting that, although along the radial coor- 
dinate direction the equation of state is above the value —1, hence being 
compatible with the equation of state for a standard scalar field, along the 
angular variables directions, depending on the heliocentric distance r^, the 
equation of state is either above, either below the value —1, hence being ei- 
ther compatible with the equation of state for a standard scalar field, either 
compatible with the equation of state for a Phantom field. 

In addition let us recall that the original motivation for considering the 
existence of either cold dark matter, either gravitational interaction cor- 
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rections is the observed flattening of galaxies rotation curves [23] which is 
equivalent to a significant correction to the gravitational acceleration to- 
wards the centre of the galaxies growing with the radial distance to the 
centre core of the galaxies. With respect to our fits we note that, for helio- 
centric distances above 70 AU up to ~ 10 pc the gravitational acceleration 
correction for the fits amax and a+ is exclusively towards the Sun decreasing 
monotonically in absolute value. As for the fit a_ the gravitational accel- 
eration correction is exclusively outwards the Sun up to the heliocentric 
distance Rmax = 95A77 AU and for the fit amm it is outwards the Sun up 
to 91 AU turning towards the Sun above this distance having a pronounced 
peak at 113 AU. Hence, if the ELA metric is interpreted as a covariant 
parameterization of cold dark matter, a extrapolation to galaxy scales with 
higher values for the central mass, tends to single out the fit o-min as the 
only fit that simultaneously allows for a positive contribution to ^cdm and 
a gravitational acceleration correction towards the central mass at large ra- 
dial distances, excluding the fit a_ and still allowing both the fits Umax and 
a+ although not justifying the contribution for VLcdm- We will develop such 
study in another work |23j . 
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A Appendix 

The coefficient values for the power series fit p!7j) obtained by a least square 
minimization method to a_ are 



02.-1 


= +1.31715 X 10^0 m 


, "2.0 


= +5.37413 X 10^ 








02.1 


= -0.00612127 


, "2.2 


= +1.07833 X 10" 
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m~ 
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"2.3 


= +1.5718 X 10-30 


, "2.4 


= -1.07371 X 10- 


43 
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"2.5 


= +3.98678 X 10"^'^ 
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= +1.59514 X 10" 
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= -4.3892 x 10-83^-7 


, "2.8 


= +3.39347 X 10- 
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-8 


"2.9 


= -9.79789 X 10-^1° m"^ 













(18) 
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and to amin are 



"2.-1 


= -1.33106 X 10^0 m 


a2.o 


= 5.69474 X 10^ , 


Q!2.1 


= -0.00573066 


"2.2 


= +2.66023 X 10-1^ m-2 


"2.3 


= 3.59912 X 10-31 mr'^ 


"2.4 


= -1.69481 X lO-^'^m-'^ 


"2.5 


= +1.18619 X 10-^'^m-'5 


"2.6 


= 9.43547 X 10"™ m,-^ , 


a2.7 


= -7.82415 X 10-^3 m-^ 


02.8 


= -6.14959 x 10-9^771-^ 


"2.9 


= 4.64438 X 10-i°9 m'^ 
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